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A New General Method of Interpolation. 

By Emory McClintock. 

A general method of interpolation will be described below, which is 
apparently both easier to prove and easier to work than the general methods 
now in use. By general method, I mean a process applicable to all cases, 
including those in which the given values of the variable do not form an 
arithmetical progression. Before stating it, it will be well to recall briefly 
the two methods hitherto prevailing. 

The first method, as we may call it for the sake of distinction, is known 

as that of Lagrange, though the credit of it has also been claimed for Euler. 

If we suppose that q>x u <px 2 , <px 3 , . . . <px k are given, where x x , x 2 , . . . x k are 

any numbers, Lagrange's formula for determining tyx, where x is any other 

number, is 

<pX = XrfXi + x 2 <px 2 + x&x a + . . . + XtfXt , 
where 



_ (x — x 2 ) (x — x 3 ). 


. . (x — X,,) 


" 1 ~~ (Xi — X 2 ) (Xi — x s ) . 

___ (x-Xi)(x—x 3 ) . . 


. . (»! — X k ) ' 

, . (x — x k ) 


[x 2 — Xi) (x 2 x s ) . , 


■ ' \X% x lc ) 



_ ( x — x ^jx — Xz) . . . (x -**_i)_ 
* _ (x k —x l ) (x k — %)... (x k — X k _i) ' 

The law by which the coefficients are formed is obvious, and the proof of the 
formula is not difficult. 

The second method, that of Newton, is as follows, using the same nota- 
tion as before. Let 

WX* — (PX, 

$& = t-^ — *U , 

%2 — *^l 

(pXo — (CX 2 

**' = -x—xV ' 

and so on, the general form being 

<px„ — <px n _i 



<piX n — 

x n x n 
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The most convenient way of recording these divided differences is to set 
down, first, a column containing the numerical values of x\, x 2 , &c ; next, a 
column containing, opposite these numbers respectively, the values of ty^, 
<px 2 , &c. ; and thirdly, a column containing ty t x 2 opposite x 2 , (piX 3 opposite x 3 , 
and so on. Adjoining this column of differences, may be recorded a second 
series of differences, <p 2 x 3 , <f> 2 # 4 , &c, calculated by the formula 

then a third series, 






and so on, the general rule being 

, fiA f m X n - 1 

9m + \X n — _ — • 

•"n ^n — m — l 

Thus, if four terms have been given, we shall have this scheme, 

X n <$>®n < Pl^« tyifin <pZ%n 



Ji\ 


Qx,. 






x 2 


$x 2 


<?>1#2 




x% 


$x 3 


<piX s 


<p 2 X 3 


"'" 4 


tyx A 


faXt 


fyXi 



To find <px, we have now to use this formula : 

<px = <px l -\-(x—x 1 ) 4>i^2 + (x— Xi)(x — x 2 ) <p 2 x 3 -\-(x — x 1 )(x—x 2 )(x — # 3 )<2>3# 4 + .... 

The complete proof of this formula is laborious.* 

Concerning the relative merits of these two methods, it may be said that 
if more than one term is to be interpolated the second is by far the best. If 
only one term is wanted, the matter is doubtful. Some writers give both 
methods; Boole, in doing so (Finite Diffs., 2d ed.), shows a preference for the 
first method, Grunert for the second. Some writers, possibly to save space, 
mention only one of the two ; Hymers gives the first, De Morgan the second. 
The first formula has often been commended for its elegance, exhibited both 
in the symmetry of its form and in the simplicity of its proof. The second, 
on the other hand, seems to be, on the whole, the easier to work in practice. 
After due examination, my conclusion is that, if all the multiplications and 



*See, for three different demonstrations, De Morgan, Calculus, pp. 550-552; Grunert, Supplement to 
Klugel's Worterbuch, vol. II, pp. 43-46; Boucharlat, Calcul, 7th ed., pp. 417-420. 
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divisions are performed by logarithms, the first method can, adopting one or 
two obvious labor-saving devices, be worked with somewhat fewer references 
to the table of logarithms than the second, but that, on the other hand, it 
requires a more frequent use of some of the logarithms when found, and 
imposes a greater strain of attention ; in short, that the second method is, for 
practical purposes, at least equal to the first. This, of course, is a matter 
which any one interested will examine and decide for himself. 

The alternative method which I have to present resembles, and may be 
regarded, perhaps, as a substitute for, the second method. The formula is 
the same, the difference lying in the manner of obtaining the quantities 
^x 2 , $ 2 x 3 , &c. I form the several orders of divided differences according to 
this formula, 



tym + \&: 



X„ — X, 



m+ 1 



In applying this formula, <£> is to be considered equivalent to <p, the given 
terms being q>x x or q>$ u <px 2 or $ x 2 , and so on. The computation of divided 
differences by this formula is simpler than by the second method, inasmuch 
as the denominator, in all the orders, corresponds to the numerator. The amount 
of labor required in both cases is the same, but the strain on the attention is 
lessened by employing a more symmetrical formula. The work done is more 
mechanical, and therefore easier. This statement may readily be tested. 

The demonstration of the new method is so simple that the process mio-ht 
be said to prove itself. Since, from the equation last given, 

we have successively, for m — 0, m = 1, &c, 

^X n = ^0^ + (X n — X\) $!#„, 
<Pl%n — fylOCt + (X n — X 2 ) $ t X n , 

q>. z x n = <p 2 x 3 + (x n — x 3 ) <p s x n , 
and so on. Hence, by successive substitution, 

^X n = 4> ^i + (X n — X t ) ^X 2 + (X n — X x ) (X n — X.,) (ptXs + (X tl — ;l\) (x a — X 2 ) 

(x n — x s ) <p s x 4 + . . . (k terms) . 
This general law is true, as proved, for all the given values of x n , namely, 
#!, x 2 , . . . x k . Assuming the same law to be true of values not given, it 
becomes a convenient formula for interpolation. This assumption is the 
same in substance with those underlying the first and second methods, and 
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the values determined in any given case by the three methods are identical. 
If we write x for x n , and <p for <|> , the formula becomes 

<px =z <px! -{- (x — x t ) ty^ + (x — Xi) (x — x 2 ) <p 2 X 3 + . . . . 

It may also be written in this form : 

<px = <px x + (x — x t ) l<p t x 2 + (x — x 2 ) {<p 2 x 3 + (x — x 3 ) [<p 3 x 4 +...]}]• 

If, for example, four terms are given, we may calculate successively 

4<3X 4 =: (x — x 3 ) <p 3 Xi, 

4> 2 x 3 =(x — x 2 ) (q> 2 x 3 + ^a) , 

4^2 —{X — Xi) ($& + 4> 2 Xs) , 

and the general formula is reduced to this, 

tyx zz tyx^ + 4 A • 
Similarly, if five terms are given, we may begin by calculating 

4 4 #5 = (X Xi) $4X5 . 

then 

4 3 T 4 —{X — X 3 ) ($3X4 + 4 4 ^g) , 

and the rest as before. The same manner of lessening the labor may be 
followed for six or more given terms. 

The results obtained by this method are neither more nor less accurate 
than those obtained by the preceding methods, being, as already stated, the 
same. In the following example, where seven values of a function of the 
sixth degree are given, the result is exactly correct. 



x x — 5 
x,— 1 1 
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x A zz 5 
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Each difference is derived from the one adjoining on the left, by subtracting 
from the latter the number at the top of the same column, and dividing the 
remainder by the corresponding difference taken from the first column. Let 
us suppose that x = 2 ; then, following the formulae last indicated, we have, 
successively, 

-^,= (2-8)^ = -— , 



720 



■4*. = (2 -6) 

fe=(2-5) 

■to =(2 -4) 

^ 8 =(2-1) 

^=(2-0) 
q>x = 5 — 5 = 



1 

240 



120/ ~~ 20 ' 
^ 20/ _ 20 ' 
' 10 W ~~ 2 ' 
^ 2 J — 2 ' 

-*+t)=-«. 



JVcfe on ^6 Demonstration of certain Formulae for Interpolation. 

When the given values of the variable are in arithmetical progression, 
interpolation is best performed by one of two well known formulae, the choice 
between the two depending upon the circumstances of the case in hand. 
These formulae are, 

<?> (x -f k) = <px + kA<px + — k (k — 1) A 2 <px + . . . , 

Li 

<p(x + k)=<px + k[A<p (x±y) + T ( A=f1 ) A 'H + F *[ a8 * O^t) 

+ ]-(^2)AV] + .... 

The notation here used is that employed in my " Essay on the Calculus of 
Enlargement" (Amer. Journal, II, 101), namely, 

k™=(k + l)k{k — l), 

k )B( = (k + 2) (k + 1) k (k — 1) {k — 2) , 



and so on ; also, 



A = E — 1, 

A = E* — E' 



-i 
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where e is the symbol of Enlargement, e a being an operation such that 
E h <px =. q> {x + h), whence 

A$x = <p (x + 1) — q>x , 

A$x = <p (x + — ) — $ (a? — - 2 ) , 

A*<px = A<l>(x + y) - A<t> (x — -i) , 

and so on. The class of differences denoted by A n $x are called central differ- 
ences. Of the two formulae above stated, the first is a simple and obvious 
case of the Factorial Theorem given in the Essay referred to, p. 143. The 
second may be derived from the same theorem, as follows. 

From the factorial theorem we have directly 

4>e = $1 + ?>e0 . A + \ <?>e0< 2 >. A 2 + ^ <?>e0< 3 >. A 3 + . . . , ( A ) 

where x 1 ® = xx , x (S) = x (x + — ) (x — y) , x (i) = x (x + 1) x (x — 1), x r "> 

— x y + y) (^ + y) (^ — "2") v ~~ t) ' and so on ' the S eneral form bein 8" 

a*"> = a? (a? + -i- m - l) (a? + -*- m- 2) . . . (x- ~ m + lV 
Also, in terms of mean central differences, 

A 2 _l- 

2.3 



4>e = 4>1 . i + <?>e0 . iA + -i- <?>e0> 2( . iA 2 + -1- 4>e0 )3 (. iA 3 + ( B ) 



Here, 



I = i-(E*+E-i), 



iA"cpx = | (A> [* + i-] + A> [or - -I]) ; 
also, 

« w =(*+t— f)(*+t— t)-(— i-+T> 

Hence, 

A rr A0 . iA + -J- A0 )2( . IA 2 + ~ A0 )3( . iA 3 + 

Since Aa' )m( = ma,' )m-1( , as was shown in the Essay referred to, 

A = iA + 0« . I A 2 + -i- )2( . I A 3 + ~- O^i A 4 + 
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But )1( =. , )8( = , and so for all odd exponents ; so that 
A = iA + y 0« iA 8 +2^0« iA 6 +. . . , 
and multiplying this by A* -1 , 

A »_ lA » + 1 )2( >lA » + 2 + J )4 <.lA W+4 + 

It follows that any even or odd central difference, or any expression com- 
posed of even or odd central differences, can be expressed, respectively, in 
even or odd mean central differences. Let 

fcE = *1 + y <?>E(P. A 2 + 2^ <?>e0«. A 4 + . . . , 
<|> 2 E = <pE0 . A + s-s <?>E0 (3) . A 3 + • • • , 

so that, observing equation (a), c^e + <£ 2 E z= <£>e. Then 

^e = a i + a 2 iA 2 + « 4 iA 4 + • • • » 
$ 2 E = «xiA + « 3 iA 8 + . . . , 

where a , a 2 , &c, a 1( a s , &c, are undetermined coefficients, and 

<?>E zr « i + a y iA + « 2 iA 2 + . . . . 
But, by (b), 

^b = ^1 . i + <?>e0 . iA + — $>e0 )2( . iA 2 + 

This determines the values of the coefficients, and it follows that the even or 
odd terms of (a) may be replaced, respectively, by the even or odd terms of 
(b). We derive in this way the following novel and important theorems in 
the Calculus of Enlargement, expressing <|>E in alternate central and mean 
central differences : 

$e = <?)1 + ?e0 . iA + \ 4>e0 (2) . A 2 + ^ ^e0«. iA 3 + 2T ^ <?>e0 ;4 >. A 4 + . . . , 
$e = &1 . i + 4>e0 . A + 4- <?>e0 )2c . iA 2 + „A; <?>e(P. A 3 + — i-j <£e(F. iA 4 + 

These general theorems can, among other uses, be applied directly to the 
matter now in hand, interpolation. Let $E =z e*, and let <tyx be the subject of 
operation ; then 

^(x + k)=^x + k iA}x + 4^ # 8) AV* + J^ k^iA'^x + . . . , 

A A • o 

^(x + k) = i^x + kA^x + ~ £ )2( iA 2 ^+ ~ ¥ Z) A^x + . . . , 

A A . O 
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These formulae are usually ascribed to Stirling, who published them in 1730. 
I have, however, seen the first and more important of the two quoted from 
Newton, 1711. They remained undemonstrated, it is said, until the discovery 
of the method of generating functions. The proof by that method, which 
may be found in Lacroix, is tedious and indirect, yet it is the only one I have 
met with, though other indirect demonstrations are certainly possible.* 

Since I = E-** =F — A, 



ix. 



iAtya? = A V (x ± -|-) =F y A n + 1 ^a 

Making this substitution in the first formula just demonstrated, and observ- 
ing that ¥ n) — k.k )n ~^, we derive the duplex formula which we set out to 
prove, 

*The difficulty does not lie in exhibiting the first few coefficients, but in proving their law. The first of 
these formulae is said by Boole to have been discussed for Smith's Prize in 1860. 



